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Motivation: Sign problem, Silver Blaze problem 
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Lefschetz-thimble path integral 


Finite-density quantum chromodynamics (QCD) 


QCD 
Fundamental theory for quarks and gluons 


Neutron star 


e Cold and dense nuclear matter 


© 2m, neutron star (2010) 


Neutron star merger 
e Gravitational-wave observations (2016~) (image from NASA) 


Reliable theoretical approach to equation of state must be developed! 


Zaco(T. 1) = [ DA Det(P(A, pa) +m) exp C Svu(). 


quark gluon 


Sign problem: Det(D(A, uq) 4- m) Z 0 at uq z 0. 
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Sign problem in finite-density QCD 


QCD & |QCD| 


ZQCD = [ra (det DP) e SYM | Z|QCD| = [ea [det P| eM, 


If these two were sufficiently similar, we have no practical problems. 
However, it was observed in lattice QCD simulation that at T = 0 


(e.g., Barbour et. al. (PRD 56 (1998) 7063)) 


nB QCD 


E > 
~My [2 -mn/3 Pa 
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Baryon Silver Blaze problem 


The curious incident of the dog in the night-time (Holmes, Silver Blaze). 
Problem: It is of great importance for finite-density lattice QCD 
to understand ng = 0 for fg < my/3. (Cohen, PRL 91 (2003) 222001) 


Current situation: For jig < m; /2, the problem is almost solved: 
Quark det. is (Adams, PRD 70 (2004) 045002) 


Det(P(A, ug) +m) | 
Det(D(A,0) +m) — 


I o60,-24. 


Re(AA) «uq 


and ess-min4(Re(A4)) = Mx/2 (Nagata et. al., PTEP 2012 014103). 


For fg > m, /2, no one knows how to understand this. 
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Method: Path integral on Lefschetz thimbles 
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Lefschetz-thimble path integral 


Sign problem of path integrals 


Consider the path integral: 


ja J «eec sta 


e S[x] is real = No sign problem. Monte Carlo works. 
e S|[x] is complex = Sign problem appears! 


If S[x] € C, eom S'[x| = 0 may have no real solutions z(t) € R. 


Idea: Complexify x(t) € C! 
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Lefschetz thimble for Airy integral 
Airy integral is given as 


3 
Ai(a) — EL (= + 2 


Complexify the integration variable: z = x + iy. 


4 


Integrand on R, and on A 
(a — 1) 
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Rewrite the Airy integral 
There exists two Lefschetz thimbles J, (c = 1,2) for the Airy 


integral: 
Ai(a) Son / sae am 
= 5 — — az]. 
= J, 2T p 3 


ng: intersection number of the steepest ascent contour Ke and R. 
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Lefschetz decomposition formula 


Oscillatory integrals with many variables can be evaluated using the 
“steepest descent" cycles 7,: (classical eom S’(z,) = 0) 


/ d”gz e- 560 = 3s R) | d”z e80, 


o 


Jz are called Lefschetz thimbles, and Im[S] is constant on it: 


ma =}, dz(t) _ eS 


t——0o dt Oz? 


I= 20) 


(Ko, R): intersection numbers of duals K, and R” 
(Ke = {2(0)|2(c0) = #2}. 

[Witten, arXiv:1001.2933, 1009.6032] 

[Christoforetti et al. (PRD(2012)), Fujii et al. (JHEP 1310), etc.] 
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Analysis: Semi-classical analysis of the one-site Hubbard model 
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One-site Fermi Hubbard model 
One-site Hubbard model: 


A = Uff, — u(fu + f). 


Fock state gives the number density immediately: 


"d 


(YT, Hidaka, Hayata, 1509.07146)(cf. Monte Carlo with 1-thimble approx. gives a wrong result: 


Fujii, Kamata, Kikukawa,1509.08176, 1509.09141; Alexandru, Basar, Bedaque,1510.03258.) 
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Analysis 


Path integral for one-site model 


The path-integral expression for the one-site Hubbard model: : 


B n A 2. UD 
Z24 — | ay (14 efe)? 78/20 
2nU R mar). 


Fermion Det 


Integrand has complex phases causing the sign problem. 


ọ is an auxiliary field for the number density: 
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Sign problem and fermion determinant 


One-site Hubbard model: 


U i 
Det lo — (» + 3 + 3] = (1 E e &CUA- Weir)” l 


Quark determinant in QCD: 


Det [yu (P 4 +m) - u] = N(A) [ [(129796:77792) (12-796), 


£j20 
where the spectrum of y4(P , + m) is 
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Silver Blaze problem for u < —U/2, u < m4/2 


One-site Hubbard model: As BU >> 1 and —U/2 — u > 0, 


U ; 
Det lo = (» + 7 + e) = (1 + EURER)" ~]. 


The sign problem almost disappears, so that J, ~ IR. 
Finite-density QCD: As 6 — oo and u < m,/2, 


Det ud, +m) — ul 
Det l4 D 4 F m)| 


The sign problem disappears by the reweighting method. 
=> Lefschetz thimbles ~ Original integration regions 


> 1 
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Analysis 


Flows at u/U < —0.5 (and u/U > 1/5) 


0.75 


0.25 


0 
Re(z) 


Figure: Flow at u/U = —1. Ix 


-0.25 


-0.5 


cR. 


N dz e- 909, 


Number density: n, = 0 for u/U < —0.5, n, = 2 for u/U > 1.5. 


(YT, Hidaka, Hayata, 1509.07146) 


Z- 
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Silver Blaze problem for u > —U/2, u > m,/2 


One-site Hubbard model: At each real config., the magnitude is 
exponentially large: 


Det lo — (» + 5 + e) —0 (e8 U+#/2)) 


This large contributions must be canceled exactly in order for n = 0. 
Finite density QCD: The situation is almost the same, since 


Det (P(A, uj) tm) _ 
Det(D(A,0) +m) — 


I er - As; 


Re(AA) «uq 


but ng = 0 for u, S my/3. 
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Analysis 


Flows at —0.5 < u/U < 1.5 
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Figure: Flow at u/U = 0 
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giu p 1. 


This value is far away from n = Im (z)/U=0, 1, or 2. 


m 


Complex saddle points lie on Im(z,,)/U 
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Curious incident of n in one-site Hubbard model 


We have a big difference bet. the exact result and naive expectation: 


95 0 0.5 1 15 ajo 


This is similar to what happens for QCD and |QCD|. 
p/U = -0.5 & py = Mr /2. 


Yuya Tanizaki (RIKEN, RBRC) Lefschetz-thimble path integral Jun 8, 2016 @ Keio 20 / 35 


Complex classical solutions 


If GU > 1, the classical sol. 
for —0.5 < u/U < 1.5 
are labeled by m € Z: 


u[U x 1/2 


U 
Ea cR By +2rmT. 


At these solutions, the classical actions become 


Sox -£ (B +3)’, 
Re (Sin — So) ~ 2m m?, 
Im S. c 2mm (5 4- 1). 


Yuya Tanizaki (RIKEN, RBRC) Lefschetz-thimble path integral Jun 8, 2016 @ Keio 21/35 


Semiclassical partition function 


Using complex classical solutions Zm, let us calculate 


oo 
Za t= > p nm. 


M=—-09 


This expression is valid for —1/2 S u/U S 3/2. 
This is calculable using the elliptic theta function: 


= 1 
-S [149 9 Heo + \a—27?m?/BU 
e ( + ) cos 27m (5 + 2 e 


Zel = 
m=1 
= H 1 2r? 
= ep |r| E+), e72 }. 
U 2 
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Number density & Lee-Yang zeros 
Lee-Yang zeros of Za: 


e o% 
. è 
r . Re 
° e 


Semiclassical study gives the correct transition! 


2 (1< WU «32. 
InZa>$ 1 (0<p/U <1), 
0 (—1/2 « u/U « 0). 


(YT, Hidaka, Hayata, arXiv:1509.07146[hep-th]) 
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Important interference among multiple thimbles 


Let us consider a "phase-quenched" multi-thimble approximation: 


Za = $e | = e 9 45(0, e 0727). 


e Lee-Yang zeros cannot appear at u/U — 0, 1. 
e One-thimble, or "phase-quenched" , result: n ~ u/U + 1/2. 


Consequence 


In order to describe the step functions, we need interference of 
complex phases among different Lefschetz thimbles. 


(cf. Particle Productions: Dumulu, Dunne, PRL 104 250402) 
(cf. Hidden Topological Angles: Behtash, Sulejmanpasic, Schafer, Ünsal, PRL 115 041601) 
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Numerical results 


Results for GU = 30: (1, 3, 5-thimble approx.: Jo, Jo U Ji, and Jo U J41 U Ja ) 
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Necessary number of Lefschetz thimbles ~ BU / (27). 


(YT, Hidaka, Hayata, NJP 18 (2016) 033002, arXiv:1509.07146[hep-th]) 
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Bonus 


Bonus: Complex Langevin method 
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Complex Langevin method 


Complex Langevin has been regarded as a sign-problem solver via 
stochastic quantization (Klauder, PRA 29, 2036 (1984), Parisi, PLB 131, 393 (1983)): 


dz(0) _ 8S 
0: Stochastic time, 7: Random force (n(0)n(0')),, = 20(0 — 0"). 
Ito calculus shows that 


x (Olen(0)) = A(O" (z,(8))), — (0 (24(8)) S (27(0)))n: 


If the I.h.s becomes zero as 0 — oo, this is the Dyson-Schwinger eq. 
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Bonus 


Complex Langevin and Lefschetz thimbles 


For any solutions of the DS eq, 


(0) = 2 374. | dee oq), 


c o 


ind, € C. To reproduce physics, d, = (Ko, R) € Z. 


So far, we ONLY assume the convergence of the complex Langevin 
method. 
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Bonus 


Semiclassical limit 
Let us take h < 1 for computing 


(Oln) = ds [ dzeS/*O(z). 


| have NO idea how to compute the LHS. However, positivity of the 
probability density and its localization around z's imply that 


"ce. 3 (Ole) S 2.508). 


RHS is 


[ dze SOM"O(z) > S Ole). 
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Bonus 


Semiclassical inconsistency 
In the semiclassical analysis, one obtains (for dominant saddle points) 


m NNNM Sieh 
Z \ Be.) l 


Co 


Co > 0. And, d, = (Ks, R) to get physics. = Inconsistent! 
(Hayata, Hidaka, YT, arXiv:1511.02437[hep-lat]) 
We show that the complex Langevin is wrong if 
e There exist several dominantly contributing saddle points. 
e Those saddle points have different complex phases. 


Open question: Recently, our knowledge on the CL method is 
developing (Aarts, Seiler, and Stamatesc, 2009; Nishimura, Shimasaki, 2015, etc.). 
Can one understand our result from this viewpoint? 
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Bonus 


Proposal for modification 
Assume that 


_ Ko, R) 2nh e- S62. 


S" (x. 


Because of the localization of probability distribution P, it would be 
given as 


P= GPs supp(P,)Nsupp(P;) = 0. 


E 


Assumption means “CL = phase quenched multi-thimble approx.” : 


(Ko, R) 


2h 
T e SGo)/h 


Sz) De): 


(0(2)) 2», 7 


o 
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Bonus 


Proposal for modification (conti.) 


If so, defining the phase function 


lS" (ze)| —iIm S(z;,)/h = 
S" (Ze) E id Xsupp(Ps) (2, Z), 


we can compute 


($(2,, Za)O(25))n 
((25, 24))q l 


This new one is now consistent within the semiclassical analysis. 
(Hayata, Hidaka, YT, arXiv:1511.02437[hep-lat]) 

Caution: Our proposal evades inconsistency, but is not necessarily 
correct. Can we improve the proposal? 


Cej 
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Bonus 


Complex Langevin study of one-site Hubbard model 


BU —30 u/U=0 


o 


o 
bo 


8[Im(z) — (u + U/2)] 
S 


= 
io 


=S  -10  -5 0 5 10 15 
BRe(z) 
Figure: CL distribution at u/U = 0 


It looks quite similar to Lefschetz thimbles. 
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Bonus 


Complex Langevin simulation 
One-site Fermi Hubbard model: 


=== exact 
2.0} HEH w/o phase 


‚9: w/ phase 

= 1.5 
o 
Ei — — — | |$*X95Emut5. .. 

0.5 

0.0 

-0.5 0.0 0.5 1.0 Ls 
p/U 


(Hayata, Hidaka, YT, arXiv:1511.02437[hep-lat]) 


Consequence 


Modified complex Langevin is not perfect yet, but it seems to point a 
correct way. 
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Summary & Conclusion 


Summary and Conclusion 


e Picard-Lefschetz theory gives a suitable framework for 
saddle-point analysis even if S(@) takes complex values. 


e One-site Hubbard model is a nice toy model to play with the 
sign problem. 


e Destructive and constructive interference of complex phases 
among Lefschetz thimbles play a pivotal role for the baryon 
Silver Blaze. 
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